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We propose optimal estimators for bispectra from excited states. Two common properties of such 
bispectra are the enhancement in the collinear limit, and the prediction of oscillating features. We 
review the physics behind excited states and some of the choices made in the literature. We show 
that the enfolded template is a good template in the collinear limit, but does poorly elsewhere, 
establishing a strong case for an improved estimator. Although the detailed scale dependence of 
the bispectra differs depending on various assumptions, generally the predicted bispectra are either 
effectively 1 or 2-dimensional and a simple Fourier basis suffices for accurate reconstruction. For an 
optimal CMB data analysis, combining all n-point functions, the choice for the excited state needs 
to be the same when computing power spectrum, bispectrum and higher order correlation functions. 
This has not always been the case, which could lead to wrong conclusions. We calculate the bis¬ 
pectrum for different choices previously discussed for the power spectrum, setting up a consistent 
framework to search for evidence of excited states in the CMB data. 


I. INTRODUCTION 

The mapping of the Cosmic Microwave Background (CMB) through the Planck [T] and WMAP [5] satellites has 
given us an incredible insight into the evolution of the Universe. It has shown us that a simple model, counting only 
six parameters, can describe the observations. Two out of six of those parameters, the amplitude and first derivative of 
the primordial power spectrum, are necessary to describe the statistical properties of the primordial fluctuations that 
source the fluctuations in the CMB. The CMB power spectrum (at least in temperature) has been measured all the 
way up to Silk damping scales, with a precision limited only by cosmic variance. In other words, with the Planck data 
in combination with small scale measurements by ground based experiments such as ACT and SPT mg, we have 
almost reached the point were no further measure of the CMB temperature will improve our current understanding 
of the early Universe. Improvement can (still) be made with the measurement of polarization, both in curl free (E) 
and divergence free (B) modes and we expect that these will indeed lower the existing constraints on As and rig in 
the coming years (as well as constraints on primordial gravitational waves). 

Theoretically, the power of the CMB lies in the fact that it potentially explores (energy) scales way beyond scales 
that can be reached with accelerators. It is therefore of great value and importance to squeeze every bit of information 
out of the CMB measurements. In particular, any deviation from a Gaussian primordial density field would provide 
a wealth of information on the early Universe as non-Gaussianity would for example tell us if the mechanism that 
produced the primordial fluctuations were driven by a weakly coupled or a strongly coupled field [g or if multiple helds 
played a role in this process. Deviations from primordial Gaussianity have already been tightly constrained; current 
constraints estimate the deviations to be at most of the order of 1 in 1000 for most well motivated shape functions. 
The limited number of modes in the GMB prevents us from lowering current bounds on non-Gaussianity much further 
and we would have to use other measures, such as the large scale structure (LSS) of the Universe, to determine the 
nature of the field [g. That being said, although profound questions about the nature of the helds will probably 
not be answered using CMB observations, a measurement of any n-point correlation function (with n > 2) would be 
immensely valuable, as these measures directly probe the interactions of the held. From an EFT perspective, the most 
valuable bispectra are the local and equilateral shapes. However, a model of the early universe with additional degrees 
of freedom can easily produce other shapes. A wide class of models produces oscillating bispectra (and higher order 
correlation functions) that break scale invariance, besides the small spectral tilt. In a recent paper we explained 
that it is hard to search for these bispectra as rapid oscillations result in a strong mismatch with the more common 
shapes, such as local and equilateral, which are only weakly scale dependent. A framework has been proposed and 
applied to search for any shape in the CMB temperature and polarization using a modal expansion, were any shape 
can be approximated using a finite number of terms on an orthogonal basis [iHia- Although this approach is very 
powerful, the used 3D polynomial basis has proven to be impractical in capturing rapid oscillations. Instead, we 
proposed a more case specific ID basis which captures such bispectra with relatively few terms. As a result, we have 
shown that in principle we can reach frequencies for the resonant bispectrum [7] way beyond the validity of the EFT 
[l4] and at the same level as the power spectrum [l5H2g and our estimator can be adopted to possible drifting [21] of 
the oscillations. 

In this paper we extend this formalism to include additional shapes. The resonant model is very specific, in the 
sense that its oscillations are logarithmic and depend only on the sum of the three momenta that connect the triangle 
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(kt = ki k 2 k^)■ The literature contains a very interesting but poorly constrained class of bispectra from excited 
states during inflation [22ll29) . Unlike the power spectrum, the bispectra are a measure of the interactions. The 
interaction Hamiltonian can be derived from EFT arguments and induces a spectrum of shapes |30| . Under the 
assumption that the vacuum is the Bunch Davies (BD) vacuum, these are all the possible shapes. In case that the 
vacuum is not the BD vacuum (for example on general grounds one might expect the vacuum state to be a mixed 
state in a multiverse scenario [29l ISU |32] ), the bispectrum is altered in two ways. First of all, the vacuum expectation 
value of the interacting fields and the interaction Hamiltonian can pick up a correction, were the correction for a pure 
state can be described by a Bogolyubov rotation. Secondly, the time integral of the interaction is altered as a surface 
needs to be specified on which the vacuum state is modified. This freedom allows for a rich phenomenology of the 
possible bispectra. Such bispectra from excited states typically have a momentum dependence of form Kj = kt — 2kj 
and include oscillations, making them an interesting target for our modified modal expansion. 

This paper is organized as follows. First, in section|^we will review the various choices of excited states that exist 
in the literature. In particular, we will try to clarify some key points regarding the decay of these signatures. We then 



II. CHOICE OF EXCITED STATE AND POWER SPECTRUM 


A. Mode equations 


The evolution of primordial curvature fluctuations C is governed by the Mukhanov Sasaki equation [34l [35] 


<+(fc2 



( 1 ) 


with the Mukhanov variable v{t) = z(t)C,(t) where = 2a^e and where primes indicate derivatives with respect to 
conformal time r. This is mathematically equivalent to an harmonic oscillator with a time dependent frequency 
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The scale factor a(T) describes the isotropic homogenous background spacetime, for example for de Sitter space 
a = and thus 

The general solution to Eq. 0 is of form 


Vk = akUkir) + l3kul{T). (3) 

The modes must be normalized so that = —i and |ofep + = 1. For the simple case of a de Sitter 

background, the positive frequency mode functions are 


Uk 



( 4 ) 


After quantization one finds that the power spectrum of the Mukhanov variable v is given by the square of the mode 
equations, i.e. Py(k, r) = |r)fe(r)p. The form of the positive and negative frequency modes Uk and is specified by the 
form of the background spacetime under consideration. The complex numbers Uk, Pk, called Bogolyubov parameters 
can be Ic-dependent. They are determined by selecting appropriate boundary conditions, and thus the state of the 
inflaton. We summarize choices that have been proposed in the literature. 


B. Review of excited states 

Before we review some of the choices found in the literature on excited states, we would like to make a general 
comment on what should be considered an inflationary initial state. It is sometimes argued that the effect of an 
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excited state will redshift away. As we move further away (in time) from the excited state, corrections should be 
suppressed; only if modifications are relatively late w.r.t. the end of inflation should we be able to observe any effects 
in the spectrum of perturbations. This is true if the modified state is applied a fixed physical time tini; indeed as 
it is a fixed time, not all modes will go through the same history and one introduces explicit scale dependence with 
any effects redshifted away from tini. This state presents a thermal bath w.r.t. to the instantaneous vacuum and 
inflation should dissipate these particles as the Universe expands. In general, the longer inflation lasts, the smaller 
the effect; if inflation lasted longer than ~70 E-folds, we expect the signal to be suppressed. These excited states 
can be considered initial states, since they can he considered as a relic of some pre-infiationary state, such as the false 
vacuum state. 

However, another possibility exists which is the consequence of considering a modification at a fixed physical scale. 
A mode crosses this scale at some point in time (which is different for each mode) and the mode is excited; later on it 
crosses another fixed physical scale, i.e. H, after which (for single field inflation) it is frozen. It is clear that each mode 
experiences the same history, besides slow-roll corrections. This is the premise of trans-Planckian physics or 

the so-called new physics hyper surface (NPH) [40] scenario. We really should not refer to such a modification as an 
initial state modification, since in some sense we are truly adding dynamics; the physical scale excites the Universe 
at every moment in time during inflation. Unlike the previous case we are not exciting each mode, but only modes 
with k < a{t)A, with A the scale of new physics. Instead, we should consider such a modification as an attempt to 
parametrize our ignorance above some physical scale, very similar to EFT. Literature suggests there is not a unique 
form of a rotated pure state in such a scenario, and several slightly different proposals have been considered. While 
each of those share the feature that they are derived from a pure state with a physical cutoff at some fixed scale, the 
details differ. 

In both scenarios, it is necessary to make sure that the energy of the particles does not spoil slow-roll inflation 
due to their gravitational back reaction on space-time. One can derive the back reaction constraints on a pure state 
rotation and one finds that the Bogolyubov parameter /3 describing the negative frequency modes, should fall off faster 
than I/fc^ |23l|4T]. These constraints are derived when summing over all modes (UV) and one way to avoid breaking 
of slow-roll inflation via back reaction is to make sure no modes are excited above some cutoff scale A. It must be 
noted that such a ‘sharp’ cutoff could potentially lead to some of the effects that we compute; if one smooths this 
cutoff some of the features that we compute could disappear. 

While UV modes need to fall off faster than I/fc^, this need not be the case for all energies. An example of this 
is provided by boundary effective field theory (BEET). In this scenario the modes are affected at some fixed time 
(the boundary hypersurface), explicitly breaking scale invariance. However, very similar to the NPH scenario, on 
that boundary there exists a (high energy) physical scale. Consequently, modes close to this energy scale are affected 
more than IR modes. The effect is thus the opposite of what one would intuitively expect for a fixed time scenario; 
the modes that are most affected are at the high energy end (and thus on smaller scales in the CMB). Such a model 
requires a UV momentum cutoff, which may be taken to be fcmax = a(Tini)A. Of course, if inflation lasts long enough 
so that physical scales in the CMB originate above this UV cutoff, no effect of the BEET state will be observable. 


1. Asymptotic Minkowski vacuum (Bunch Davies vacuum) 

The standard choice of vacuum is defined as the lowest energy asymptotic initial state (i.e. without particle 
excitations). This choice is known as the Bunch Davies (BD) vacuum. In the infinite past, r —>■ oo and every mode is 
deep within the horizon and much smaller than the curvature scale. One can then use the Minkowski initial conditions 

lim Ufc(r) = (5) 

V2fc 

which selects the positive frequency modes, i.e., ak = 1, bk = 0. In other words, the BD vacuum is the equivalent of 
the Minkowski vacuum in de Sitter space. 


2. Constant Bogolyubov rotation 

The algebraically simplest possible modification of the Bunch Davies vacuum is a ^-independent Bogolyubov ro¬ 
tation, i.e. Cfc = Cl, bk = C 2 with constant complex coefficients so that |afcp -I- = I. As we shall review below, 

such a rotation only changes the amplitude of the power spectrum, and is therefore indistinguishable from the power 
spectrum in Bunch Davies. However, as we shall later show, the momentum dependence of the bispectrum is modified. 
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Physically this is because a non-zero Bogolyubov rotation implies that the modes are not in their Bunch Davies vac¬ 
uum state, and therefore contain particles with respect to this state. The bispectrum measures interactions between 
these particles, and can therefore be enhanced. 


3. a-vacuum at a fixed scale and NPH scenario 


An common inflationary state proposal (see e.g. [38) 1 is to excite each Fourier mode when its wave length becomes 
equal to a new fundamental scale A (for example the Planck scale). The model would thereby avoid to describe scales 
outside the range of validity of the theory. As argued earlier, such proposal is really dynamical and therefore we 
should not refer to it as an initial state modification. It assumes, the initial conditions for each mode are set at a 
^-dependent time Ti^iik) such that k = a(Tini)A. We set these initial conditions at a length scale and time scale so 
that the mode is deep within the horizon (this is possible since A i7), so that the mode equations are to good 
approximation of Minkowski type and we have 




(Ani) 




( 6 ) 


A natural assumption is that high energy effects should be suppressed by powers oi H/A. The Bogolyubov param¬ 
eters can then be put in the perturbative form [42] 


111 , i^(A„i(fc)) , 

\ak\ = 1+1/- \-0 

H{Ti^i{k)) 

\Pk\ = X - ^ - VO 


A ^) ’ 
A^ 


( 7 ) 


Our guiding principle here is to construct initial conditions that are exactly scale invariant in de Sitter space. 
Therefore we assumed that the magnitude of the Bogolyubov parameters depends on k only implicitly due to H 
(when breaking de Sitter invariance by making H a function of k), and x and y are k-independent real constants, 
naturally of order 1. In addition, the relative phase of the positive and negative frequency solution must be chosen 
scale invariantly, since it is not an overall phase of the quantum mechanical state and therefore observable in the 
power spectrum, as we shall review in the next section. In de Sitter space we have Tini = The unique scale 

invariant vacuum choice at this order in H/A is therefore given by the Bogolyubov parameters 

afc = (l + (8) 


Here we chose the phases of the modes so that the a-modes start at Tini in phase 0 and the /3-modes start in the 
fc-independent phase cj) (when inserting Eq . ([^ in Eq. (©)• 

To compare our result with that of Ref. |40p we rewrite the mode equation in the form v = N{k){u}~ -\- b{k)u\) and 
obtain 

b{k) = (9) 

in agreement with the Bogolyubov parameter of the new physics hyper surface (NPH) scenario. They also note that 
in this model the exponential factor appears to be necessary to avoid non-localities to order H (while here it arises 
due to scale invariance). Both the NPH scenario and the trans-Planckian argument are based on the same principle; 
there is some physical scale (NPH: A, trans-Planckian: Mp) beyond which we assume each mode rotates off the Bunch 
Davies vacuum, i.e. is excited. Since this is a fixed scale, very similar to the Hubble scale during inflation, every mode 
will experience the same history, up to slow-roll corrections due to the broken de Sitter symmetry. Any mode that 
exited the horizon before set scale is not affected and is assumed to be in the vacuum state; however these modes are 
unobservable as they are beyond our current horizon. 


4- Fixed time vacua and the BEFT scenario 

A true initial state modification sets the initial conditions at a fc-independent time Tini, for example the earliest 
time in which we can trust general relativity. The initial conditions are then explicit functions of k, breaking scale 
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invariance. In this scenario, corrections to the Minkowski vacuum should come as powers of of the physical wave 
number = fc/aini divided by a cutoff scale A. The Bogolyubov parameters are then 


\ak\ = 1 + 2/ 
\Pk\ = X 


k 


OiniA 

k 

^ini A 


o 

-o 


OiniA / ’ 

k 


^ini A 


We assume that all modes start at Tjni at a fixed k-independent phase. Therefore we find 


afe = 1 + 2/ 


k 


niA 




I3k = [x 


niA 


„i(fcTi„i + 0) 


( 10 ) 


( 11 ) 


Unlike in the NPH case, oscillations in the power spectrum here also arise in pure de Sitter space, where H is constant. 
The Bogolyubov rotation in so-called boundary effective field theory (BEFT) [43] was first parametrized in Ref. [40] 


b{k) = X 




( 12 ) 


Compared with the previous case, the non-BD vacuum here explicitly breaks scale dependence. Since there is both 
a fixed time and a fixed scale, different modes, frozen when they cross the horizon, have not experienced the same 
history. As we discussed above, the BEFT scenario can generate large effects on UV scales, which is a result of the 
assumption that on the initial hypersurface, modes that have a momentum closer to the set UV cutoff scale, will be 
more affected. Because of the form of the Bogolyubov rotation, we expect this modes to be best constraints using 
smaller scales, for example those probed by LSS at high redshift. However, given the back reaction constraints, modes 
k > at some time Oini are not excited. Although not explicit in this equation, we should keep in mind, if this 

is too early, the comoving cutoff will eventually redshift outside of our current horizon and we will not be able to 
observe these effects at all. 


C. Shape of the primordial power spectrum 

As discussed by Ref. [40| the Bogolyubov parameter b{k) = |&(fc)|e“^^) defined in Vk = N{k){uk + b{k)u%) immedi¬ 
ately leads to a power spectrum of the curvature perturbations of form 

P{k) ~ Pbd (1 + 2|5(fc)| cos(a(fc) -I- ijj)). (13) 

A constant Bogolyubov rotation thus only results in a renormalization of the power spectrum (at lowest order in slow 
roll). For the NPH initial conditions discussed before we get 

Pbd (l + +'*)) ■ (14) 

where H(T\^i(k)) oc We thus obtain logarithmic oscillations in k. The BEFT scenario gives 

pBEFTQ^) ^ Pbd (\j^2x-^cos{-^ -(15) 

where Hint is constant at Tjni, so we obtain linear oscillations in k. The phase tp is arbitrary and fc-independent. 

Generally, if the solution to the e.o.m. for the inflation degree of freedom can be written as a super position of 
positive and negative frequency modes, and if there is scale dependent phase, oscillations in the power spectrum will 
appear 131101111]. Consequently, if the breaking of scale invariance is at most slow-roll, models are hard to distinguish 
at the level of the power spectrum. Higher order correlation spectra proof a way out, since they are sensitive to the 
dynamics of the field and as such, allows us to discriminate between models. 
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III. SHAPE OF THE PRIMORDIAL BISPECTRUM 


A. General remarks for arbitrary Hint 


1. Leading order bispectrum 


The bispectrum at tree level for a Gaussian initial state is given by |3S] 

(CkiCksCka) = / ^^([CkiCk2Ck3.^fl('r)]) 

= -27^e J dri(CkiCk2Ck3^i(T)), (16) 

with Hi the interaction Hamiltonian. The Hamiltonian is theory specific as are the solution to the equations of motion 
for C (which by itself are derived from the varying the action w.r.t. to C). We start with general remarks about this 
calculation, before evaluating the bispectrum for the simplest cases. 

Without further details on Hi we can consider solutions for C that are rotated w.r.t. to the usual Bunch Davies 
vacuum, i.e. are in an excited state. Each rotated (^, will now contain some particles, and the density is determined by 
the Bogolyubov parameter /3. Unlike the power spectrum, a constant /3 already changes the momentum dependence of 
the bispectrum (while for the power spectrum a fc-independent rotation would lead to a renormalization of the power 
spectrum). This is sourced by the fact that interactions take place during inflation and the strength of this effect is 
determined by the interaction Hamiltonian as well as the total time that inflation lasts. More precisely, the presence 
of particles and the interaction affect both the shape and amplitude of the primordial curvature 3-point function. 

If the solutions of the e.o.m. for C, are a superposition of positive and negative frequency modes as explained in 
section [II A[ then given the Hamiltonian to cubic order in C, and 3 <C’s for the 3 vertexes we then generally find that 
the 3 point function can be written as (to lowest order in the /3 and with a ~ 1) 


(Ck 3 Ck 2 Ck 3 ) 0(/3°)+ci(fc)(5(^k,) J dr 


■E- 

j 


AKjT 


-I- c.c. 


0 (/ 3 "), 


(17) 


with Kj = kt — 2kj and the amplitude of ci depends on the strength of the interaction and its scale dependence is 
such that for the chosen interaction (and thus the power a") the overall scaling of the bispectrum is as 1/fc®^. Note 
that since this is a connected 3 point function, that the only surviving contraction are those with 3 distinct (ingoing) 
momenta, the other 3 momenta serve as conservation of total momentum. 

Note that for n = 0 we have ci(fc) oc 1/k^. A more complicated interaction Hamiltonian will generally lead to 
different power of a, and hence different pre-factors in powers of Kj. This is crucial, as different powers of a also 
dictate the maximal enhancement [231124] ; for n different from 0 one has to partially integrate to perform this integral 
and therefore higher power of Kj will appear in the denominator. The higher the power the more powers it can 
‘cancel’ in the numerator, e.g. 


while 


lim 


sin ujKj 


= w, 


lim 


1 — cos iuKj 

K'f 




(18) 


(19) 


etc. Since w 3> 1 the enhancement in the collinear limit is dictated by the power of a. 


2. Choice of initial time surface and k dependence of the Bogolyubov parameter 

To evaluate the shape via Eq. 0), there are two choices to be made that were not present in the power spectrum 
calculation, and that influence the resulting shape. These are the initial time surface Tini(fc) and the momentum 
dependence of the Bogolyubov parameter (which can now depend on either k^ or kt for example). 


If the cutoff time is scale independent, one explicitly breaks scale invariance and this argument no longer holds. 
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We first discuss the choice of initial time. It is clear that the above integral diverges for finite momenta if we take 
the conformal time integral all the way to negative infinity (except for the part inside 0(/3*^)). From back reaction 
constraints it was obtained Ham] that /3 should fall off faster than 1/fc^. Here that is achieved by making sure that 
no modes are excited for k > aA Unlike the power spectrum, we now have to be specific about our initial time. We 
shall see that the choice here is crucial and determines the scale dependence of the bispectrum. We consider several 
possible choices, in line with our discussion of the effects on the power spectrum. The most obvious choice for the 
initial time is to place a finite time cutoff tq, independent of the momentum k. However, the following argument 
suggests another, momentum dependent initial time surface. For all positive frequency modes, we should choose 
an initial condition that returns the bispectrum of single field slow-roll model. The contribution from the positive 
frequency modes goes like 


(CkiCk 2 Ck 3)/30 ~ C2(fc)(5(^ k.,) J dr + c.c.] . ( 20 ) 

For the Bunch Davies choice, one typically gets terms that go like l/kt- However, when we apply a modification at 
some time Tini we get additional scale dependence introduced by the cutoff; these corrections are of the form 


sin/ctTini 

h ’ ^ ^ 

and a possible choice of Tjni could be m: jkt x A/i7, which makes sure this term vanishes and all positive frequency 
modes lead to the solution associated with BD initial conditions. One can argue for a third choice. As we have 
seen above, once we add negative frequency solutions, we introduce a problem since an excited vacuum would give 
infinite contribution to the 3-point function at early times. From this perspective it is the negative frequency mode 
that should determine the initial time surface, so one should choose a cutoff that renders this direction a ‘special’ 
direction and Ti = 1/fci x AjH. This is slightly different from the choice made earlier {Tini{kt)) and would generate 
an excited state earlier since kj < kt, so its effect should be bigger. We would like to point out that with all three 
choices, Tjni = const., Tini(fct), rini(/ci), the computation relies on a sharp cutoff. Such a cutoff might not be completely 
natural. In case the boundary conditions are softer, one would expect that some of the features computed here might 
be damped. 

The momentum dependence of the Bogolyubov parameter is also crucial. Even when using a constant fc-independent 
Bogolyubov parameter, one gets oscillations in the bispectrum. This is in contrast to the situation in the power 
spectrum. For example, the characteristic logarithmic oscillations in the NPH power spectrum explicitly require 
H = H{k), i.e. H is broken by slow-roll and the associated Bogolyubov parameter obtains a scale dependent phase. 
In previous works it was generally assumed that the Bogolyubov rotation was constant and hence did not lead to 
additional scale dependence; eventually, one would like to compare the observations of the power spectrum that 
of higher order correlations functions (e.g. the bispectrum) and for that reason we will also consider a Bogolyubov 
parameter that is /c-dependent, consistent with the choices made for the power spectrum. If the Bogolyubov parameter 
is explicitly momentum dependent, as is the case in BEFT, the natural momentum dependence is l3{ki), where ki 
is the momentum of the corresponding leg in the excited state. If on the other hand depends on the momentum 
only implicitly hy H — H{k), as in the NPH scenario, we use the total momentum of the three-point interaction, i.e. 
P{H{kt)). The time cutoff Tini in the NPH scenario is naturally chosen with the same momentum dependence. To 
summarize, for the NPH scenario we will use Tini(fct) and P{H{kt)), and for the fixed time BEFT [43] scenario we use 
Tini = const, and P{ki). This choice also guarantees that the resulting shapes are at most effectively two-dimensional in 
their modal expansion to be proposed below, independent of the powers of k coming from the interaction Hamiltonian. 


B. Single field slow-roll non Bunch Davies shapes 


We now specialize our above discussion to the single field slow-roll interaction Hamiltonian |5S] 


Hi 


d^x 





( 22 ) 


which captures the dominant non-Gaussianity contribution in a single vertex. While single field slow-roll is well known 
to give small non-Gaussianities, we nevertheless get insights into the possible oscillating fc-dependencies induced by 
non-Bunch Davies vacua. We are using de Sitter mode functions, as in the case of the power spectrum. While 
corrections to the mode functions arise when iJ 7 ^ 0, as in the case of the NPH scenario, such corrections are slow-roll 



for this interaction Hamiltonian 


suppressed. From Eq. (|20|), to linear order in /3fe, the bispectrum correction due to the non-Bunch Davies vacuum is 

1,2 jL 2 1,2 
/vi«,2fi.3 iji 


^ 2 tjQ fO 


k- 


-hc.c.. 


(23) 


As reviewed above, there are different possible choices for the fc-dependence of Tini and f3k, which we will discuss next. 


1. Constant Bogolyubov parameter, constant or k-dependent tini 

The algebraically simplest result is obtained for a constant Bogolyubov parameter /? = and a constant initial 
time Tini. In this case, the integral evaluates to 




-(2^)3<53(^ 


dije klklkl 


E 


COs((/>) 


— cos{ujKj -I- (/>) 

-rfK, 


(24) 


where the oscillation frequency is formally given by the initial time, uj = tq. In the case of rini = const., in principle 
also the term with three legs in the BD state is modified, and acquires an oscillation in cos{ujkt). This form of 
oscillatory shape was discussed in Ref. [33] where the CMB multipole bispectrum was calculated and an optimal 
estimator was proposed. 

If the initial time is chosen as rini oc the shape is modified to be 


^(C/:iCfc,Cfc3)nBD2 = -{2TTfS^{Yl 


4ij6 kfk^ki 


E 


COs{(j)) 


cos{{uj^) + (j)) 


(25) 


which changes the directional dependence of the oscillation in fc-space. 


2. NPH scenario 


Now we assume an NPH vacuum, i.e. /3fc = (x^) e (dropping a possible constant phase (j) for notational 
simplicity) and Tini = A/{Hkt). We find 


A/r A A \ 1 klk'^kl fH\ cos{jj) cos( ^^ ( J -f 1)) 

A(C£.‘.£.f4>.BD3 =-(2») « (L‘-))(^n(2tJ)^(»A 


k^K, 


, which has the correct limiting behavior for Kj —^ 0, i.e.. 




kik2k3 


fcf 


(26) 


(27) 


So far we have not taken into account the breaking of de Sitter symmetry in the background, that is necessary to 
obtain the logarithmic oscillations in the NPH power spectrum. In the case of the power spectrum, this could be taken 
into account by making /? implicitly depending on k via H. To be completely consistent, one should calculate the 
NPH bispectrum in the slow-roll expansion, to take into account both the slow-roll corrections to the mode functions 
and the slow-roll corrections in the vacuum rotation Pkt ■ Here we only take into account the latter of these effects 
and obtain 


, , k ,o klklkl^f cos(j^logfct)-cos( 7 j^(^ + l)logfct) 


In the collinear limit we then obtain 


Jim^A(Ck4Ck2Ck3) oc 5(^fo)^^5] 

3 


k^Kj log kt 
sin(A/iJ*log kt) 




fc3 


(28) 


(29) 


which resembles resonance non-Gaussianities. 
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FIG. 1: The shape of Eq. (311 on the diagonal = £2 = ^3 for 4 different frequencies. The constant determines the overall 
scaling with i, and superimposed oscillations appear on top of the BAO. 



FIG. 2: The shape of Eq. (311 in the collinear limit £3 = 500 = £1 -\- £ 2 - Unlike the diagonal there are no visible oscillations, 
even when you increase the frequency. 


3. Fixed time vacua with fik 


In this case we explicitly break scale invariance in the initial conditions by considering (3^ as in the BEFT scenario 
P5] . For a fixed time vacuum at Tini with Bogolyubov parameter j3kj = we get linear oscillations of 

form 


A(Cs,Cs,C^3)nBD5 = 


klklkl 

MpV2 n(2fcf) 



cos (Tinifcj) 


-COS (Ti„i(fcj +Kj)) 

k^K, 


(30) 


As is the case in the power spectrum, we find that for fixed time initial states the non-BD contribution of the 
bispectrum is more suppressed the smaller the physical momentum k/aini is compared to the energy scale A (and 
therefore the larger the scale on the CMB sky). Note that in principle there could be additional contributions from 
the boundary itself as discussed in [15] . 


C. Exact calculation of example CMB multipole bispectra 


In this section we evaluate exact CMB multipole bispectra for some of the shapes discussed in the previous section. 
We use these exact calculations to obtain an intuition about what can be seen in the CMB as well as to verify our 
modal expansion proposed in the next section. 

The simplest case is the shape derived in Eq. (24) with corresponding normalized bispectrum (assuming here phase 
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0 = 0 ) 


, , , A:?fcp| 1 - coswiiT, 

^ 1^2 ^3 j 


kp<, 


( 31 ) 


We can use the original re-parametrization proposed by Fergusson and Shellard [S]: 


ki = ka = fc(l — /3) 

k 2 
ks 


kb = -fc(l + a + P), 


kc = — a + /3). 


(32) 


Note that ki -\- k 2 + k^ = kt = {a + b + c)k = 2k, hence k = kt/2. The parameters have the following domains 
O<A:<oo,O<0<l and —(1 — /3) < a < (1 — /3). The volume element can be computed through the determinant 
of the Jacobean, i.e. Det Jij = k^, and dkidk 2 dk^ = k^dkdadjd. Unlike the resonant shape, the example shape above 
now explicitly introduces an non-factorizable form in the a-jd plane. This is not a real problem, since the original 
integrals over the Bessel and transfer function already require sampling both in a and in fd in one loop. One can 
actually propose a very similar parameterization to avoid this, but because there is really no computational advantage, 
we decided to stick to the original parameterization. It is not hard to show that the shape is then given by 


S{q,a,ld) = 


1 (1 — cos 2 g(l — a)a;) 


a^bc 


(1-a) 


perm. 


(33) 


In the shape above, a can be I and numerically this can lead to some issues which are easily resolved by explicitly 
writing down the a, 6 , c —>■ I limit. As it should, that limit is finite and does not diverge; It is enhanced in the squeezed 
limit but not beyond local m- 

In Fig. we show the projected shape one the equilateral axis = £2 = ^3 divided by the constant bispectrum 


^constant 


{2£i + 1) (2^2 + 1) ( 2.^3 + 1) _it + 3 £t_ 


(34) 


From the figure, it appears the shape is dominated by the constant, which simply leads to BAOs in the bispectrum. 
As we increase the frequency of the oscillations, the shape roughly shows the same BAO structure, with superimposed 
oscillations, very similar to what happens to the power spectrum. Also note that for w = 50 partial cancellation seems 
to occur between the oscillating part and the BAO, resulting in overall lower amplitude. 

In Fig. we show the collinear limit, with £3 = 500 and £±+£2 = £ 3 ,- This is the same limit in which the primordial 
spectrum maximizes, i.e. corresponding to Kj —>■ 0. Again we should not be concerned with absolute value of the 
amplitude, as it is arbitrary. Unlike the equilateral limit there are no clear visible oscillations on top of a smooth 
spectrum. Even when you increase the the frequency, no clear features appear, suggesting that the collinear limit is 
indeed dominated by a constant, i.e. oj. In the collinear limit, the peaks appear when either one out of 3 £ gets large, 
i.e. £2 "C -^3 ~ £1 (and symmetric equivalent) and as the frequency is increased £2 ^ £1 = £ 3 / 2 . 

We will consider another example, in case the Bogolyubov rotation and the initial time Ti„i both depend on kt as 
in Eq. (28 1 we have^ 


I ^cosa;logfct-cosa;(I-^)logfct 
B[ki,k2,k3) - ^ 


klklklh^ fc2(l-^)logfct 


We find 


I coswlog 2g — cos(a;(l — 2a) log 2 ( 7 ) 
b[q,a,jd) = -~™—+ perm.. 


a^bc 


(I - a) log 29 


(35) 


(36) 


with the limiting behavior 


lim S(q,a,jd) = lo sin(w log 2q) + perm. 
a,b,c^o a^bc 


( 37 ) 


^ Inside the log, kt should be dimensionless, so we choose some pivot scale k* = 1 Mpc 
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FIG. 3: The shape of Eq. ( |35| l on the diagonal = £2 = (-i for 4 different frequencies. Compared to the shape in Fig. 
the constant has now been replaced with a cosine oscillating in logfet; this leads to a shape that resembles something like the 
resonant bispectrum in the equilateral limit. 



FIG. 4: The shape of Eq. ( |35[ ) in the collinear limit ^3 = 500 = Small features appear on all scales, with strong features 

in the squeezed limit. 


We show the result for several frequencies in the equilateral limit in Fig. The dominant contribution is now set 
by the cosine oscillating in logfcj. It therefore resembles the resonant bispectrum in this limit, however the shape is 
modified since the smooth scale dependence now is a function of Kj. 

In Fig. I^we show the collinear limit for = ^l+i 2 = 500. For the frequencies considered the bispectrum changes 
sign, but hardly changes amplitude; it does however appear to contain rapid oscillations in the squeezed limit. 


IV. MODE DECOMPOSITION AND ESTIMATOR FOR THE CMB BISPECTRUM 


A. The case for modal expansion 


A template shape, called the enfolded shape, has been proposed to recover the collinear behavior of the shapes 
above [2l]. It was shown that as we increase the frequency of the oscillation, the overlap decreases between the 
template and these non-BD bispectra. Also the proposed template does not have the correct scaling in the squeezed 
limit, making it less useful for for example bias in large scale structure. 

The CMB multipole correlator between two bispectra B and B' is given by 


Corr(i3, B') 


_ 1 _ 

norm(B)norm(S') ^ Ci^Ci^Ci^ 


(38) 


We first consider the collinear limit, before calculating the full (all tj correlator below. To show that the enfolded 
template is in fact useful in the collinear limit, we have computed the overlap in the collinear limit between the 
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FIG. 5: The enfolded template as proposed in [ 21 ] compared to the shape of Eq. (31 1 and Eq. (35l. It has strong visnal 
resemblance and a more quantitative comparison reveals that these shapes overlap 97 — 99% in this limit. It also shows that 
although the enfolded captures the collinear limit, it does worse in the sqneezed limit, where the visual resemblance is less 
apparent. However, since it is finite in this limit, it does not ruin the total overlap. 


enfolded template and Eq. (31) on the collinear axis above (i.e. £3 = 500 = £1 + £ 2 )- We find that the overlap on the 
collinear axis increases from 80% for w = 50 to 99% for ui = 10000. The reason why the overlap goes up is that as we 
increase the frequency the collinear enhancement becomes much more pronounced. In Fig. we show the collinear 
limit of the enfolded and non-BD shape of Eq. (311 (with some arbitrary normalization). It is clear that the features 
in this limit are very similar, specifically the enhancement in the limit £1 = £ 2 . It is also clear that this template does 
not do well in the squeezed limit (i.e. when £1 <C £ 2 , £ 3 ), which is relevant for using the template for the computations 
of the non-Gaussian bias in large scale structure. We find similar overlap for the shape in Eq. (35), with overlap 
almost independent of frequency ^ 95 — 97%. 

To show that there is a need for an improvement over just a simple enfolded template, we computed the overlap of 
the full shape up to £max = 2000 with the enfolded template for various frequencies, both in primordial k space and 
in CMB multipole space. We use the primordial shape correlator [111 [48] 


F{Si,S 2) = j dVSi{ki,k2,kz)S2{ki,k2,kz)W{ki,k2,k-i). 


(39) 


and {B, B') = F{Si, S2) / {F{S{)F{S2)) and the weight proposed in |9| >V(fci, fc2, ^ 3 ) = k^^ . We show both the 

priordial and the CMB correlator between the enfolded template and Eq. (|3H as a function of w and for two values of 
the phase in Fig. First of all, we find that as we increase the frequency the overlap decreases as was shown earlier 
Ref. [24]. Secondly, the enfolded shape is most adapted to a phase <() = 0. Any deviation from that generally leads 
to a worse overlap, as was pointed out in Ref. [25j : since the phase is a free parameter, we would severely limit the 
extend of parameter space we are sensitive to if we would rely on the enfolded template only. Thirdly, after projection 
the overlap improves as oscillations are washed out. This was also found in Ref. 

It is worth pointing out that the primordial correlator is a function of fcmax- Even though our fcmax is set to 0.7 Mpc“^ 
for the CMB multipole bispectra, in order to obtain somewhat reliable overlap in primordial space, the primordial 
correlator has to be cutoff around the Silk damping scale fcmax fcsiik = 0.1 Mpc“^. Depending on their functional 
form, some oscillating shapes are less sensitive to the value of fcmax- One should be aware of this difference when 
comparing the results presented here with those found in Ref. [24] ; there the shape considered has a scale dependent 
cutoff. Consequently, the overlap did not depend on /cmax and the effective frequencies are somewhat higher and 
therefore the overlap smaller. In addition, for the more complicated shape of Eq. (35), there is a fc-dependent phase. 
As a result of the scale dependent phase of this shape we find a very small overlap (< a few percent) for most 
frequencies and phases. We only investigated the overlap of this shape in primordial space, but although some of the 
oscillations will be washed out by projection, we do not expect large improvements of the correlator in CMB multipole 
space. 

Our aim is to squeeze every bit of information out of the CMB and therefore we prefer a representation of the 
theoretical shape that is > 99%. For that purposes we require a modal expansion, which will be the topic of the next 
section. 
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FIG. 6: Primordial (3D) and CMB multipole (2D) bispectrum correlation of the enfolded template and the exact non-BD shape 
Eq. (311 as a function of frequency cj for 2 values of the phase. Although the overlap in the collinear limit matches well with 
the enfolded template, rapid oscillations off the collinear limit result in a poor overlap. After projection some of the oscillations 
wash out and the overlap improves. We also investigated the overlap with the second shape (Eq. (|35[)), and found much worse 
correlation due to the scale dependent phase, establishing a strong case for a more accurate reconstruction. 


B. Mode decomposition 

As we we have seen above, non BD bispectrum shapes are generally of the form B{kt, Kj), with appropriate 
prefactors of single momenta ki that are not an obstacle for separability. Our goal is to find efficient expansions of 
our primordial shapes in terms of separable basis functions, exploiting the symmetries of the primordial shapes. This 
will allow us to construct efficient estimators and make signal-to-noise forecasts in the following sections. 


1. Effectively 1-dimensional shapes 


The simplest shape of interest is Eq. (311 which is only a function of Kj. In this case, we can use a 1-dimensional 
expansion similar to that in Ref. [7]. In general any function f{x) can be expanded in a Fourier series 0 ^ 6 ™^, 
where x in the present case is some function of the momenta ki. The resulting expansion will be useful for our purpose 
if the exponentials are separable, i.e. if a; is a linear function of momenta ki and not a more complicated function of 
these. Therefore we can expand Eq. (31) as 


3 N 

b{Ki,K2,K3) = J2 E ' 

j—1 n— — N 


with expansion coefficients 


Cn — 


AKd 


rK'7 




B{K,)e 


dK,. 


(40) 


(41) 


The range of integration over the 1-dimensional variable Kj corresponds to the physically allowed values of Kj, 
i.e. 0 < Kj < 2kmax, where kmax is the maximal single momentum for which we want to describe the primordial 
bispectrum. The bispectrum must be periodic on the space of integration, which we achieve by multiplying with a 
generalized gaussian window function (see Ref. H). Fig. [7] shows an example expansion of the primordial bispectrum 
in Eq. © with (w = 800, (/) = 0). We see that our expansion in terms of oscillating basis functions is well suited to 
represent both the fast oscillations and the enhancement for Kj ^ 0. 


The CMB multipole bispectrum for the shape in Eq. (31) is given by 


N 

= dxx^ Cn{A';j:{x)A'^{x)B'i^^{x)+A‘^{x)B‘\{x)A'^{x) + B‘\{x)A'^{x)A‘^{x)), (42) 

_Ar 


'i^-N 


where 


Ai{x) = ^ J dkk‘^An{k)ji{kx)Ai{k), Bf^{x) = ^ J dkk'^Bnik)je{kx)Ai{k), (43) 

























14 




FIG. 7: Primordial shape f{Kj) 


cos{ujKj ) 


(top) and corresponding modal coefficients c„ = a„ + (bottom) for to = 800. 


and An{k) = and Bn(k) = Fig, [s] shows the bispectrum calculated with our expansion using 600 modes 

compared to the exact result of section [III Cf 


2. Effectively 2-dimensional shapes 


For shapes B{kt,Kj) like Eq. (25) we construct an efficient 2-dimensional modal expansion. The general 2- 
dimensional Fourier expansion, of a function f{x,y) where x and y are two independent variables, is 


f{x,y) = 

mn 

We can therefore expand the bispectrum as 

Bikt,K,,K2,K3)=Y,Y. 


Cm.n^ ^ 6 


. 27Tnkt 

AiC, Afc* 


j mn 

The Fourier coefficients are given by the 2-dimensional integral 


1 


f'K" 


AKjAkt JJ^min 


B{Kj,kt)e * e~"^^dKjdkt, 


(44) 


(45) 


(46) 


where 0 < Kj < 2kmax and Skmin < kt < Sk^ax- Defining m' = and n' = and expressing the bispectrum as 
a function of single momenta ki, thereby exploiting the relation of kt and Kj, we find 

B(ki,k2,k3) = E Cmn [-^m'+n' — m'-\-n' (fca) -|- Ajn'+n'{kl)B 

— m'+n' {k2)A,^'+n'{k3) (47) 

mn 

+ B — m'-\-n' -^m'-\-n{^2)-^m'-\-n' (^s)] ; 

i27rafc i27rafc 

where Aa{k) = ^ and Ba{k) = —. Here s,t are the k scalings that have to be adjusted to the exact shape 

under consideration. 

The CMB multipole bispectrum is then given by 


bhhh = J dxx'^ c™„Hj;^,+„,(a;)Hj^,+„,(a::)HL^„,^_„,(a;) -b (perm, -b m -m). 


(48) 
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FIG. 8: Comparison of the exact bispectrum results from section [III C| and the expansion proposed here, using 600 modes. Top: 
Equilateral limit = I 2 = (- 3 ■ Bottom: Collinear limit I 3 = 500 = d-i + (. 2 ■ The results coincide as expected. 


The time critical part of the CMB bispectrum is the calculation of the functions A\^{x) and By choosing 

the size of the integration domains AKj and Akt identical or a multiple of each other, one can reduce the number 
of functions that have to be calculated considerably. While our 2-dimensional expansion is computationally more 
demanding than a 1-dimensional one, it is still a considerable improvement with respect to a general 3-dimensional 
modal expansion for arbitrary shapes. In addition, the use of oscillating basis functions makes it optimally suitable 
for the oscillating shapes created by non-Bunch Davies initial conditions. 


C. KSW type modal estimator 

The KSW estimator for the 1-dimensional expansion in Kj is similar to the one for kt in Ref. [7]. It is given by 
f = ^ (Scuh + 'S'lin) where, using complex quantities, the cubic term is 

^ 2 ^^ f dVLAl^(r,fi)B_n(r,n), (49) 

n 

and the linear correction (to take into account noise and partial sky coverage) is 

S'^“ = -3^c„ J dr y dD [R„(r, h) (^^(r, n))-b 2^„(r, h) (R„(r, n)A„(r, h))], (50) 


where (...) are the usual Monte Carlo averages over simulated gaussian maps. For the shape in equation Eq. (31) the 
filtered maps are given by 


A„{r,n)= ^a)imAiir)Ytmin), B„{r,n) = ^{C ^a)emB^(r)Yim{n), 


(51) 


im 


with Ai{r), Bi{r) given in Eq. ( [4^ . The norm N of the estimator will be calculated in the next section. The 
expansion factors Cn must be evaluated for a sufficiently dense sampling in oj and (j). 
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For the 2-dimensional expansion, the estimator is given by 


:;*cub 


E< 


r'^dr 


—m'-\-n' (j"•) ; 


(52) 


and 


= J r'^dr j dVt[B_rn'+n'{r,h) + 2Arn'+n'{r,n) {B_rn'+n'ir,h)Arn'+n'{r,n))]. (53) 

mn 

The computationally most demanding part of the KSW estimator, the computation of the filtered maps, does not 
change fundamentally when going from the Id to the 2d expansion. Only the multiplication of these real space maps 
in the estimator scales quadratically with the number of modes. However this part is subdominant in the calculation, 
making the estimator computationally feasible. 


D. Fisher forecast 


In this section we calculate the Fisher matrix, which provides the signal-to-noise ratio as well as the normalization 
of the KSW estimator, for the simple 1-dimensional shape in Eq. (31). For a CMB experiment with noise power 


spectrum Ni and sky coverage fraction /s/cy, the Fisher matrix for bispectra indexed by i,j (here discriminating 
between frequencies , ujj ) is 


Bij — fsky 


I 1 I 2 I 3 






(54) 


where the reduced CMB multipole bispectrum is here given by eq. 42 Ci = Cf^^ + Ni is the power spectrum. 


and is the usual geometric factor. We normalize the bispectrum in analogy with other oscillating shapes in the 

literature as 


B^{ki,k2, fca) = 


A|/nl 

{kik2k^f 


1 — COS ojKj 


k]K, 


(55) 


and we assume A| = 9.04 x 10 The diagonal Fisher matrix value Fa as well as the corresponding standard 
deviation c/ = of optimal measurement with the estimator of the previous section is plotted in figure(top). 
For low w, the cosine is nearly constant and therefore the bispectrum Eq. (55) goes to zero. From the Fisher matrix 

which is plotted in Fig. (bottor 


one can also calculate the correlation coefficient corr(/i,/j) = 


Dm) for some 


yF~F~, ^ _ 

test frequencies. It gives an impression of the necessary sampling in frequency space when doing the estimation. As 
one might expect from the visualization of spectra in previous sections, there is a large overlap in wide ranges of 
frequencies. 

For the resonant bispectrum [7] cr increased as a function of ui. For the resonant bispectrum, the increasing error is 
caused by the rapid oscillations. For the bispectrum considered we also expect rapid oscillations (off-collinear), but 
the produced bispectra from excited states are naturally enhanced in the collinear limit with and enhancement linear 
in w for the bispectrum of Eq. (31). If the scaling would hold after projection, we would expect that a would drop 
accordingly, so as l/w. However, as was pointed out in [24], the presence of oscillations and because of projection one 
looses some of this enhancement. It was argued through an analytic computation that one would in fact loose the full 
enhancement ui in Ref. |53|, but here and in Ref. [23], we roughly find that tr cx l/logw. In other words, the error 
decreases as you increase the number of oscillations. 

For realistic values of (H/Mpi) this shape is most likely unobservable in the CMB. An exception would only occur 
if the fundamental scale of the underlying high energy theory is similar to the energy scale of inflation, and well below 
the Planck scale, which increases uj and thus the remaining enhancement. Our main objective in this paper was to find 
suitable expansions for the typical fc-dependencies of non-BD shapes, and the present shape is merely an illustration 
of the most simple case. More complicated Hamiltonians with higher derivative operators would lead to potentially 
observable signals H5] . even for small values of oj. Although the scale dependence of these bispectra is more 

complicated, the overall scaling is preserved, including the presence of features and enhancement in the collinear limit. 
The tools presented here to reconstruct the simplest case, are straightforwardly extended for more complicated shapes 
and we will include the search for these spectra in the data in future work. 
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FIG. 9: Top: Fisher matrix value Fu (left) and corresponding standard deviation (right) as a function of frequency uj. Bottom: 
Correlation matrix corr(/i,/j) as a function of frequencies uji,ujj. 


V. CONCLUSIONS 

In this paper we developed an optimal set of estimators for bispectra from excited states during inflation. Although 
the detailed scale dependence of the bispectra can differ based on several assumptions, one can distinguish two key 
properties that are very general; first of all, the predicted bispectra are enhanced in the collinear limit, and secondly, 
the bispectra contain oscillating features. Both of these phenomenological properties make it hard to use existing 
techniques to build fast, optimal estimators. Similar to the reconstruction of the axion bispectrum in [7], we used the 
fact that all these bispectra are effectively at most 2-dimensional, allowing for a much faster reconstruction using a 
simple I or 2-dimensional Fourier basis. We tested our expansion, and showed that we can achieve good overlap with 
the exact result, using a limited number of modes, even for highly oscillating bispectra. 

We also calculated the overlap with the enfolded template, which has been used to search for bispectra from excited 
states in the CMB. We found that in the collinear limit the predicted bispectra have large overlap with the enfolded 
template. In Ref. [24] it was shown that the overlap decreases as a function of the frequency of the bispectrum (for 
an NPH pure state rotation); we find that although this is true, in the collinear limit itself, the overlap generally 
increases as a function of frequency. As was pointed out in e.g. Ref. [^, the enfolded template does not have 
the correct scaling in the local limit (making it less useful when considering effects on the halo bias), however the 
squeezed contribution does not diverge. The enfolded template has poor overlap anywhere else, in particular when 
the bispectrum has a non-vanishing phase. This establishes a strong case for the reconstruction approach presented 
in this paper. 

The ultimate experiment combines measurement of all the n-point correlation functions. Hence, in this paper 
we have stressed the need for consistency when predicting effects from excited states. Although the differences are 
typically of order slow-roll in the power spectrum, at the very least we should be aware of them. For the bispectrum 
they can lead to much larger differences and it is crucial that the correct estimator is applied when searching for the 
associated signal of a given model. 

Together with the optimal estimator for the resonant bispectrum, we have now completed a framework that allows 
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us to constrain a large class of bispectra that were difficult to constrain with existing methods. In this paper we 
focused on a very simple model, with only a single field and no non-canonical kinetic terms and a quadratic potential. 
In a more complete picture EilllT], more complicated bispectra are predicted. The framework presented here can 
easily be extended to capture those shapes as well. We hope to report on this search in future work. 
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